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ABSTRACT 



We report moderate resolution 3-5 fim spectroscopy of the nucleus of 
NGC 1068 obtained at 0.3" (20 pc) resolution with the spectrograph slit aligned 
approximately along the ionization cones of the AGN. The deconvolved FWHM 
of the nuclear continuum source in this direction is 0.3". Four coronal hues of 
widely different excitations were detected; the intensity of each peaks near ra- 
dio knot C, approximately 0.3" north of the infrared continuum peak, where the 
radio jet changes direction. Together with the broadened line profiles observed 
near that location, this suggests that shock-ionization is the dominant excitation 
mechanism of the coronal lines. The depth of the 3.4 /xm hydrocarbon absorption 
is maximum at and just south of the continuum peak, similar to the 10 /im silicate 
absorption. That and the similar and rapid variations of the optical depths of 
both features across the nucleus suggest that substantial portions of both arise in 
a dusty environment just in front of the continuum source(s). A new and tighter 
limit is set on the column density of CO. Although clumpy models of the dust 
screen might explain the shallowness of the silicate feature, the presence of the 
3.4 /im feature and the absence of CO are strongly reminiscent of Galactic dif- 
fuse cloud environments and a consistent explanation for them and the observed 
silicate feature is found if all three phenomena occur in such an environment, 
existing as close as 10 pc to the central engine. 

Subject headings: dust, extinction — galaxies: active — galaxies: nuclei — infrared: 
galaxies — galaxies: individual(NGC 1068) 
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Introduction 



The unified model of active galactic nuclei (AGNs) explains the most basic observed 
properties of AGNs in terms of toroids of dust and gas, whose orientations to our line 
of sight determine the observed characteristics of these objects. The infrared spectral 
energy distributions (SEDs) of AGNs betray the presence of dust in their nuclei. Some 
insight into the properties and distribution o f the dust has been 



and modeling of the 9.7 /im silicate feature (iMason et al. 



Roche et al. 


2006 


2007) 


Mason et al. 


2004 


2007) 



2006 



Honie & Beckert 


20C 


7; 


band ( 


Dartois et al. 


2004 



20071 ). Notable findings include differences in the spectral profile of the 



lOum silicate feature between galaxies with differing amounts of obscuration ( IRoche et al. 



20071 ). a nd evidence that t 



microns (Jaffe et al. 



Raban et al 



2004 



le toroids are rather small, only a few parsecs in diameter at 1 



Packham et al. 



2005 



Poncelet et al. 



2006 



Tristram et al. 



2007 



20091 ). At shorter wavelengths many AGNs are known to have nuclear point 
sources, whose angular extents are diffraction-limited (e.g., at ~0.1" resolution as viewed in 
the i^'-band). 

The dust clo se to an AGN is generally expected to be accompanied by large quantities 



of molecular gas (IKrolik fc Begelman 



19861 ). so it might be assumed that medium or high 
resolution infrared spectroscopy would result in the detection of molecular absorption lines 
there. Surprisingly, probably the easiest interstellar infrared molecular band to detect in 
absorption, the CO fundamental (f=l-0) band, whose most prominent lines are located 
at 4.5-5.0 yum, has not yet been fo und toward the nucleus of the prototypical and nearby 



14.4Mpc, 1^^= 72 PC 



resolution (iLutz et al. 



TuUv 



20oX and at 20 km s ^ resolution (Masor 



19881) Seyfert 2 galaxy, NG C 1068, at b oth 



et al. 



20 km s"^ 



2006r). Line emission 



from molecular hydrogen in portions of the central few parsecs ( iMiiller Sanchez et al. 
is the only direct evidence for gas phase molecules there. 



20091) 
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A clumpy torus was first suggested by iKrolik fc Begelmaru ( 19881) to account for dust 



Nenkova et al. 



fl2002h to explain the 



surviving close to AGNs and has been invoked by 
behavior of the silicate feature with orientation and the spectral breadth of the far-infrared 
le mo del is supported in some specific cases by fits to the SED of an AGN (e.g., 



emission. T 


le mo 


Honig et al. 


2008) 



20081 ). In NGC 1068 the presence of prominent dust absorption bands both 



to longer and shorter wavelengths than the CO fundamental band suggest that some CO 



absorption ought to be pr esent. In Ga 



of 0.4, the value found by 



Mason et al. 



a ctic d ense cloud material, a silicate optical depth 
(120061 ) toward the nucleus of NGC 1068, implies 



an column density of 0.7 x 10^^ cm ^ and a CO colu mn density of 1.0 x 10^*^ cm 



( Roche fc Aitken 



1984 



Bohlin et al. 



1978 



Lee et al. 



19961 ). which should be easily detectable 



in the fundamental band. Toward the diffuse cloud in front of Cyg 0B2 No. 12, which has 
Ay ~ 10 mag, the 10 fim silicate absorption is of comparable strength to that in NGC 1068 
and the 3.4 /im absorption band, generally assurn ed to arise from stretching of C-H bonds 



in solid-state hvdro carbons (|Sand: 



in NGC 1068 (e.g.. 



Whittet et al. 



brd et al 



1997 



19911). is consi derably weaker than that band 



Mason et al. 



2004J ). N evertheless several (narrow) 



lines in the CO fundamental band have optical depths of 0.3 (iGeballe et al. 



19991). 



One way for CO to be present but difficult to detect is if the CO lines are so broad 



that they blend toget 



Gebhardt et al. 



ler, so that individual 



(l2000h and 



ines cannot be seen. Using different methods 



Bian &: Gul (120071 ) have each estimated the mass of the nuclear 



black hole in NGC 1068 to be ~2 x 10^ Mq. The dimension of the outer of the torus- like 



clouds observed by 



Jaffe et al. 



Raban et all t(]Q% 



( l2004j ) and recently re-observed by 
using mid-infrared interferometry is 2 pc. Together these values imply an orbital speed of 
225 km s^^, slightly less than half the spacing of lines in the CO fundamental band. Thus, 
if absorbing CO is located at the outer edge of the torus, each CO transition would produce 
an absorption 450 km s~^ wide, which is comparable to the CO line spacing. In that case 



and assuming that the CO rotational states are thermalized to at least a few hundred 
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Kelvins, the M band spectrum would show a broad and shallow depression covering most 
of the band and would not be easily detectable in previous data sets. Any CO orbiting 
considerably farther from the nucleus, however, should be readily detectable in individual 
lines, if present with sufficient column density. For example, gas orbiting at a radius of 
10 pc would produce lines with full widths at zero intensity (FWZIs) of ~200 km s~^. 

The high signal-to-noise ratio attainable in the M band for NGC 1068, has allowed 
us to conduct a new and more sensitive search for CO absorption against the M band 
continuum at the nucleus of this galaxy. We also have obtained a 3-4 /im spectrum of 
the nucleus in very good seeing and with high enough S/N to begin to discern the spatial 
distribution of the 3.4 /im feature within 0.4" (30 pc) of the nucleus. In the Galaxy 
this band, due to the C-H stretch i n ethyl and methyl groups attached to h ydrocarbons, 



is observed only in diffuse clouds (jPendleton et al. 



1994J : 



Chiar et al. 



20021 ). Recent 



observations of the spatial behavior of the 9.7 /xm silicate feature toward AGNs have 



provided some information on the dust distribution in the centra . 



( iMason et al. 



2006 



Rhee fc Larkin 



2006 



Roche et al. 



2006 



2007 



few tens of parse cs 



Young et al. 



2003), in 



several cases reveal ing extended dusty s tr uctures perhaps assoc i ated w ith material fueling 



the central engine (jPackham et al. 



2003). 



Miiller Sanchez et al. 



(120091 ) have found stronger 



evidence for fueling of the NGC 1068 AGN, based on high angular resolution observations 
of H2 line emission within the central few arc-seconds. More detailed observations of the 
3.4 /im band and a more sensitive search for CO absorption promise a complementary view 
of the material in the center of NGC 1068. 

The 3-5 /^m region also includes a number of high-excitation forbidden emission lines. 
Such "coronal" lines in AGNs have been a subject of interest during the last few decades 
for several reasons. First, their high ionization energies (>100 eV) imply that they must be 
excited by shocks and/or hard photons from the AGN, with negligible contribution from 
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hot stars. If the possible sources of excitation can be disentangled, for example through 
observation of the locations at which various lines arise, the line strengths and ratios may 
yield valuable information on the nature of the AGN continuum in otherwise inaccessible 
regions of the UV/X-ray spectrum. Second, if they are emitted in the inner regions of 
the narrow line region (NLR), or even associated with the putative torus of the unified 
model, they can be used indirectly to study those structures. Third, the widths, profiles 
and wavelength shifts of the lines may map extreme conditions such as shocks and outflows 
in regions close to the central engine, which may be spatially resolvable in nearby AGNs. 

Optical spectroscopy from Hubble Space Telescope (HST) has provided detailed 
information on the morphology of the emission line regions and complex velocity fields 
of the ionized ga s associated with radio structures close to the AGNs (for NGC 1068 see 



Axon et al. 



19981 ). Although efforts have been made to understand the infrared spectra of 
AGNs, because of the lack of high spatial resolution information on a well-defined sample, a 
consensus on the location and origin of the infrared coronal lines has not yet been reached. 
In this respect, NGC 1068 is a useful laboratory because of the prominent coronal line 
spectrum that it displays at optical, near-infrared and mid-infrared wavelengths. Whereas 
the prominent coronal lines in the optical are largely from Fe, the 1.0-2.5 fim region includes 
several other species, notably Si, Ca and Al, and the 3.0-5.0 /^m range includes lines from 
Al, Si, Ar, and Na. 

The new 3-5 fim observations and their reduction are outlined in ^ and the strengths, 
distributions and profiles of the emission lines and molecular and solid state absorption 
features are described in §3.2l and §3.31 respectively. Comparisons with previous observations 
are mostly found in §4l The implications of these findings for our understanding of the 
environments producing the coronal lines and solid state and possible molecular absorptions 
are discussed separately in ^ 



Fig. 1. — Orientation of the 0.2" slit used to obtain the 3-5 /im spectra. N orth is up and 



east is to the left. Edges of the slit (vertical lines) are 



(119971 ). a 4" X 4" field containing the O III image from 



contours of 



Muxlow et al. 



overlayed on Fig. 1 o 



Macchetto et al. 



Capetti et al. 



fll994h and the 6 cm 



(119961 ). Radio source SI, as designated e.g. in iGallimore et al. 



(120041 ). is coincident with the near-infrared continuum peak. 
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Observations and Data Reduction 



Spectra of the nuclear region of NGC 1068 and of a nearby calibration star, HR 
691 (FOV, L=4.7), were obtained at the United Kingdom Infrared Telescope (UKIRT) 
on Mauna Kea on 2006 Au gust 10. The facility 1-5 fim imager/spectrograph UIST 



( jRamsay Howat et al. 



2004J ) was used, with its 0.2" (2-pixel wide) slit centered on the 
continuum peak and oriented north-south, as shown in Fig. [H With that slit width UIST 
delivers resolving powers, R, of 1400-2400 in various parts of the 3-5 /im region. The 
orien tation roughly corresp onds to the position angle of the ionization cones in the central 



0.4" (IGalhmore et al. 



19961 ). Rows of the array are separated by 0.10" along the slit. Table 1 



summarizes the observational parameters. 

The observations were made in superb seeing conditions and through a photometric 
and relatively dry sky (1.5 mm of precipitable H2O). Three spectral regions were covered, 
2.9-3.6 fim, 3.6-4.2 /im, and 4.4-5.3 /im. Total exposure times were six minutes in each 
of the above spectral intervals, and the telescope was nodded 12" north and back, along 
the slit in an ABBA pattern whi ch was repeated eve r y two minutes. Lo ng-slit spectra and 



imaging of H I Br 7 at 2.17 fim (iTamura et al 



1991 



Davies et al 



19981 ) have detected no 



line emission at the offset location. 

HR 691 was observed in each spectral interval just prior to NGC 1068 as a flux 
calibrator and to remove telluric absorption lines from the spectra of NGC 1068 during data 
reduction. The differences in mean airmass between the pairs of observations of NGC 1068 
and calibration star were each one percent or less. The intensity profiles of the continuum 
from HR 691, measured along the slit have full widths at half maximum (FWHMs) of 
approximately three pixels 0.3" (equivalent to 22 pc at the distance of NGC 1068) in all 
three wavelength bands observed. Assuming that the intensity distribution observed along 
the slit is circularly symmetric, we estimate that 55 ± 10 percent of the signal from HR 691 
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fell out of the slit. 

Data reduction employed Starlink Figaro routines to extract spectra in individual rows 
of the UIST detector array from coadditions of several subtracted pairs of images that had 
been fiat-fielded in the UKIRT reduction pipehne. Following manual removal of spikes, 
spectra from pairs of adjacent rows were combined to produce spectra covering various 
0.2" X 0.2" regions along the slit. These were ratioed by the extracted spectrum of the 
calibration star. Wavelength calibration, obtained from telluric lines in the spectrum of 
the calibration star is better than 0.0001 /xm. All wavelength scales in this paper are in 
vacuo and not corrected for the systemic velocity of NGC 1068. The flux calibrations of 
the reduced 2.9-3.6 /im and 3.5-4.1 /im spectral segments agreed to 3.5 percent in the 
wavelength interval of overlap, attesting to the excellent weather conditions. Much of the 
4.4-5.4 iim region that lies outside of the 4.5-5.0 /xm interval is unusable or of little value 
because of low telluric transmission and is not shown here. In addition, we do not include 
the reduced spectrum in several narrow wavelength intervals in the 4.80-5.00 //m region, 
where deep telluric H2O lines severely degrade the spectra of both the galaxy and the 
calibration star. 

In the 3-5 fim interval the spectrum of the calibration star contains several prominent 
absorption lines of hydrogen in the Brackett, Pfund, and Humphreys series. Of these only 
Br a falls in a relatively clean portion of the spectrum and can be removed by interpolation 
across the observed feature. For the other prominent H I lines we employed a Kurucz 
R=4,000 model spectrum of Vega with line strengths somewhat diluted in an attempt to 
match the H I lines of an FOV star. This technique proved successful in removing the Pf P 
and Hu 11-6 pair at 4.65-4.67 //m but did not work well for simultaneous removal of the 
Pf 7 and Pf 6 lines in the short wavelength 2.9-3.6 fim segment. Small wavelength intervals 
near those lines are omitted in the spectrum presented here. 
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Results 



3.1. Dimension of Nuclear Source 



Full widths at half maximum (FWHMs) along the slit were determined for both 
NGC 1068 and the calibration star from the coadded spectral images The results are shown 
in Table 1. In all three wavebands the nuclear continuum source in NGC 1068 had a FWHM 
of approximately four pixels along the slit (compared to three pixels for the calibration star) 
and was thus partially resolved. A simple Gaussian deconvolution, using the stellar FWHM 
as the instrumental value yields intrinsic FWHMs for NGC 1068 of 0.27"-0.29" in the NS 
direction (corresponding to a characteristic NS "radius" for the continuum emission of 
~10 pc). This appears qualitati v ely co nsistent with the adaptive optics imaging near these 



wavelengths by 



Marco fc AUoinl (1200011 as we. 



higher resolution imaging of 



Gratadour et al. 



1 as w ith the more recent and considerably 



(120061 ) ■ Because the central source is partly 



resolved in the present data, it is possible to crudely investigate the spatial distribution 
of the 3.4 /im absorption feature along the slit, in addition to accurately determining the 
distributions of the emission lines. 



3.2. Atomic Lines 



Figure [2] shows the 3-5 /im spectrum extracted from the central ten rows (i.e., a 
0.2"xl.O" region of NGC 1068). The expected locations of hydrogen recombination lines and 
fine structure lines are indicated for a heliocentric syst emic radial ve 



In addition to the broad 3.4 /im hydrocarbon feature (iBridger et al. 



1997 



Mason et al. 



emic radial ve 


ocity of 1137 km s ^ 


1 

(iBrideer et al. 


1994; 


Imanishi et al. 



2004J ). four high excitation fine structure (coronal) lines are seen in 
emission: [Alvi] at 3.660 /im, [Siix] at 3.9357 /im, [Arvi] at 4.5295 /im and [Navil] at 
4.6847 /im (rest wavelengths). Of these the [Alvi] and [Navil] lines are new detections. 
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: NGC 1068 0.2"xl.O" 
- 20060810 UKIRT/UIST 


1 1 1 1 1 1 




hydrocarbon 






CO band center 




[Ar VI] [Na VII] 

Pf(S 



4 
4.5 



4.6 4.7 4.8 4.9 

Observed Wavelength (/U.m) 



Fig. 2.— Fig. 1: 3-5 fim spectrum of the central 0.2" x 1.0" (EW x NS) region of NGC 1068. 
Wavelengths of atomic lines of interest are shown, as are the the three components of the 



3.4 /im hydrocarbon feature (jPendleton et al 



19941 ) and the band center of CO, all at the 



systemic redshift of 0.003793. The noise level varies but can be estimated at each wavelength 
from the typical point-to-point fluctuations of the continuum in the vicinity. For simplicity 
the flux scale assumes that both NGC 1068 and the calibration star are unresolved point 
sources and that slit losses for them are equal; this is approximately correct for the continuum 
from NGC 1068, but not for the much more extended line emission. 
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A tentative detection of the latter was reporte d by 



were previously found by 



Marconi et al 



fll996h and 



Lutz et al 



Lutz et al 



(2004). The other lines 



(|2000| ). The coronal lines 



have FWZls of over ~1000 km s ^ and the peak of each line is blue-shifted relative to the 
systemic velocity of NGC 1 068 by ~ 300 km s~ ^ , as has been found previously for other 



coron al lines in NGC 1068 (iMarconi et al. 



1996 



Lutz et al 



2000 



Rodriguez- Ardila et al 



20061 ). No hydrogen or helium recombination lines are detected in Fig. [2J The marginal 
emission feature near Pf 7 has a profile and a velocity centroid that are inconsistent with 
all of the other lines, leading us to doubt its reality. There is some marginally significant 
spectral structure in the vicinity of Br a, near 4.06 fim, but no obvious sign of line emission. 

The coronal line that is detected at the highest signal-to-noise ratio is the [Si ix] line at 
3.95 fim. Si^^ has an ionization energy of 303 eV, making it one of the most highly excited 
species detected in NGC 1068. The next best coronal line for study in these data is the 
[Alvi] line at 3.67 fim, which is more than three times weaker, but is situated at shorter 
wavelengths than other lines and in a region of high atmospheric transmittance so that the 
ambient background is relatively low and the noise slightly lower than at 3.95 fim. The 
ionization energy of Al+^ is 154 eV, approximately half that of Si"*"^. Figure [3] shows both 
the continuum flux density and the fluxes of these two lines, in 0.2" x 0.2" apertures along 
the north-south slit. The line fluxes peak approximately 0.3" north of the continuum peak 
and emission extends northward to 0.8". Little or no line emission is present south of the 
nucleus. The peaks of the line emission from both [NaVll] (ionization energy 172 eV) and 
[Arvi] (ionization energy 75 eV) also are offset a few tenths of an arc-second to the north 
of the continuum peak, similar to [Siix] and [Alvi], but their detailed distributions are 
less certain, because of difficulties in accurately removing telluric absorption features and 
accurately defining the continuum level in the case of the [Ar vi] lines, and because of the 
weakness of the emission in the case of the [Na Vll] line. 
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Fig. 3. — Distribution of 3.95 fim continuum, [Siix], and [Al vi] line flux in 0.2" sections and 
0.2" steps along the 0.2" wide slit. Typical ±lcr uncertainties in line strengths are shown; up- 
per limits, designated by downward arrows, are 2a. The angular resolution (FWHM = 0.30"), 
derived from the proflle of the standard star, is indicated. 
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NGC 1068 [Si IX] 3.9357^m 
UKIRT 20060810 0.2"x0.2" apertures 




1000 



LSR Velocity (km/s) 



Fig. 4. — Velocity profiles of the [Siix] line in 0.2" steps along the slit. Continuous hues 
denote spectra at the center and to the north; dashed lines indicate spectra to the south. 
Numbers in parentheses are continuum offsets. The y-axis is flux density in a 0.2" x 0.2" 
aperture and assumes an estimated sht loss factor of 2.2 for the calibration star. The systemic 
radial velocity and velocity resolution are shown. 
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Figure m shows the velocity profiles of the [Siix] line along the north-south slit. The 
blue-shifted peak seen in the composite spectrum of Fig. 1 is present at all positions from 
the center to 0.8" north. At these locations the velocity of peak line emission is shifted 
between -200 and -400 km s~^ relative to the systemic velocity of NGC 1068. To the south 
no line emission is seen except perhaps 0.4" south where a weak line centered near the 
systemic velocity may be present. At all of the locations where the line is clearly detected, 
the velocity profile is asymmetric with an extended blue wing. The wing is most prominent, 
essentially forming a shoulder, in the central and 0.2 " north spectra. Overall the line 
emission appears to extend approximately to -1000 km s~^ and +400 km s~^ from the 
systemic veloci ty. This ra n ge is similar to that seen in other coronal lines in the 1-4 /xm 



region, e.g., by 



LutzetaL 



(12000f ). 



The line fluxes within the UIST slit given in Table 2 were determined in various ways. 
Coronal line fluxes derived solely from Fig. [2] will be systematically incorrect because the 
spectrum has been fiux-calibrated using a point source, but the line emission is more 
extended along the slit than is the continuum. The spectra in Fig. H] and summed fluxes 
for the [Silx] and [Alvi] lines in the table have been corrected for slit losses, using the 
spectrum of the calibration star. For [Arvi], and [Navil], the signal-to- noise ratios are too 
low to determine fluxes by summing their values at locations along the slit. Therefore we 
have assumed similar distributions for these lines as [Silx] and [Alvi]. Then their fluxes in 
the slit can be estimated by scaling the fluxes derived from Fig. [2] by 0.60, the ratio of the 
fluxes of the [Si ix] and [Al vi] lines in Table 2 to the values for them derived from Fig. [21 
For Br a and Pf 7 upper limits were derived solely from the spectrum in Fig. [21 
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Table 1. NGC 1068 Observing Log 10 August 2006 



Wavelength Exp. Time R NGC 1068 FWHM'^ HR 691 FWHM'^ 
range seconds A/AA arcsec arcsec 

2.9-3.6 fim 360 1400 0.41±0.02 0.31±0.02 

3.6-4.2 fim 360 2400 0.42±0.02 0.30±0.03 

4.4-5.3 /im 360 2000 0.41±0.03 0.29±0.03 

'^Measured along the slit at 3.4 /im, 3.9 /im, and 4.6 fim. 



0.04 I 




0.1 

-1 -0,5 0.5 

N offset (arcsec) 



Fig. 5. — Optical depth of the 3.4 /xm hydrocarbon absorption feature (this paper) and the 
9.7 yum silicate absorption feature (iMason et al.ll2006l ) across the nucleus of NGC 1068. Error 
bars are ±la. 
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Table 2. 


, NGC 106; 


8 Nuclear Line 


r luxes 


Line loniz. Energy 


WdVcicli^Lii 


JTlLiA 111 OllL 




eV 


rest /xm 


10-16 W m-2 


[Alvi] 


153.8 


3.660 


0.45±0.10 


H I 8-5 


13.6 


3.7405 


<0.3 


[Si IX] 


303.2 


3.9357 


1.5±0.15 


H I 5-4 


13.6 


4.0522 


<0.8 


[Ar vi] 


75.0 


4.5295 


2.5±0.7 


[Navil] 


172.1 


4.6847 


1.1±0.3 



'^Fluxes of coronal lines are derived from Fig. H] and in 
the cases of [Ar vi] and [Na vil] by scaling values derived 
from Fig. [2] as described in the text. The upper limits to 
the recombination lines are derived from Fig. [21 Upper 
limits and uncertainties (in parentheses) are 2a; uncer- 
tainties are mainly due to uncertainty in the continuum 
level. 
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3.3. Solid-state and molecular absorptions 



The 3.4 fim hydrocarbon absorption feature is the only solid-state feature present 
in the 3-5 fim spectrum. It extends from 3.35 fim to perhaps 3.57 fim with a broad 
absorption maximum near 3.41 fim and a shoulder that begins at 3.47 fim. Assuming 
that the contribution of the 3.43 fim component drops to zero at 3.50 fim, the absorption 
maximum of the shoulder feature occurs close to 3.50 /im. After accounting for the redshift 
of NGC 1068, the above wavelengths and the overall prof ile are in good agreern ent with 



those observed for the 3.4 /im feature in Galactic sources (jPendleton et al 



1994h 



The optical depth of the 3.4 fim feature has been determined by fitting a linear 
continuum across the 3.10-3.35 fim interval (excluding the 3.29-3.33 fim region) and the 
3.57-3.60 fim interval. The maximum optical depth of the 3.4 fim feature shown in Fig. [2] 
is 0.072 ± 0.004. 

Using the same technique we also have measured the optical depth of the 3.4 fxm 
feature in 0.2" steps within ±0.4" of the infrared continuum peak. The results are shown in 
Fig. [5l together wi t h the measurements of the optical depth of the 9.7 fim silicate absorption 



from 



Mason et al. 



(120061 ) measured in 0.4" steps steps over a wider region. At and south of 
the continuum peak the optical depth of the 3.4 fim feature is approximately 0.08, while 
at the northern locations the depth is noticeably less, approximately 0.06. This pattern is 
similar to that seen in the silicate feature. The implications of this are discussed in §3 

The 4.5-5.0 ftm portion of the spectrum allows a new search for the presence of the 
fundamental band of carbon monoxide (CO), which should manifest itself as individual 
absorption lines if the CO line width is less than the line spacing (~500 km s~^ or as broad 
depressions on either side of the band center if the CO line profiles are comparable to or 



broader than the spacing. 



Mason et al. 



(120061 ). observing at a resolution of 20 km s ^, did 



not detect narrow CO lines. As discussed earlier CO absorption lines arising in gas orbiting 
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the AGN outside the characteristic radius of the continuum emitting region should be no 
wider than 200 km s-^ (FWZI). 

No signs of absorption by individual lines can be seen in the current spectrum, and any 
broad depression between the band center and the long wavelength cannot be deeper than 
about 1% of the continuum. Close examination of the spectrum shows that the continuum 
in the 4.8-4.9 fim portion of the spectrum, corresponding to P branch transitions with 
10<J< 20, may be depressed by roughly that amount relative to adjacent portions of 
the spectrum. This could be interpreted as evidence for absorbing CO in dense gas at 
temperatures of several hundred degrees and with very high velocity dispersion. For an 
isothermal absorbing slab, lines of CO at this temperature with FWHM's of 500 km s~^ 
and central depths of 1% would produce such a depression and would imply a CO column 
density of ~ 1 xlO^'' cm~^ and a gas particle density of order 10'' cm~'^. The corresponding 
portion of the R branch of CO is not available for comparison due to telluric absorption 
and contamination from the [Ar vi] line. Because the evidence for this broad CO absorption 
is marginal we regard the above column density as an upper limit. 



4. Comparisons and Constraints 



4.1. Atomic lines 



Both the offsets of the peaks of the [Siix] and [Alvi] lines from the peak of the 
continuum and the e xtents of the line emission are somewhat less than those reported by 



Marconi et al. 



(119961 ) for the coronal lines that they measured along the ionization cone at 
position angle 31 degrees. However, their measurements were made at roughly five times 
lower angular resolution than those here. Thus the present determination of the offset 



should be more accurate. The distributions of the other two coronal lines along the NS slit 
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are not as well known as those of [Siix] and [Alvi], but they appear to be similar. 



Thompson et al. 



(I2OOII ) imaged the 1.96 /xm coronal [Sivi] line at 0.2" resolution, 



finding it to originate mainly to the north-northeast of the the co ntinuum source. 



current measurements of the [Si ix] line are largely consistent with 



Thompson et al. 



but failed to detect the fainter line emission that their data show to the south. 



'he 



(20011 



Prieto et al. 



(I2OO5I ) imaged the [Sivil] 2.48 fim line emission at 0.10" resolution. Their figures show it 
to be mostly in a clumpy structure centered on the continuum pe ak. This appears t o diffe r 
markedly from the distributions of both the [Si vi ] line imaged by iThompson et al. (2001) 



Marconi et al. 



(119961 ) and by us. The 



and the [Silx] line observed spectroscopically by 
accuracy of the co-registration of their images is not available. In the present data the lines 
and the continuum were measured simultaneously and the result is clear. The ionization 
energy of [Sivil] lies between those of [Alvi] and [Silx]; thus it would be surprising if the 
intensity distribution of the [Sivil] line relative to that of the adjacent continuum were not 
similar to the lines observed here. 

The new measurements of these lines put significant constraints on the size of the 
high-ionization line emitting region. The detected [Si ix] fl ux in the UIST slit is about 



3.5 times smaller than that reported by iLutz et al. 



(I2OOOI ) using the Infrared Space 



Ob servatpry (ISO ) with a 14" x 20" aperture, and three times smaller than that reported 



by 



Marconi et al 



(Il996l ). who used a 4.4" slit width. However, assuming that the [Silx] 
line-emitting region through which our observations have sliced is roughly circularly 
symmetric, the total flux in the [Silx] region is ~5.3 x 10^^^ W m~^, in agreement with 
the flux measured by IS O and, within the uncertainties, with that by Marconi et al. Using 



adaptive optics imaging 



Prieto et al. 



(I2OO5I ) fund a rather clumpy distribution of the [Si vil] 



2.48 fim line emission extended over ~1" but one that does not prefer any particular 
position angle, so the assumption of circular symmetry, while likely incorrect in detail. 
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should allow a reasonable estimate of total line flux. 



Following the same line of reasoning, the [Alvi] line, which has a very similar 
flux distribution along the UIST slit as [Siix] (see Figure [3]), would have a total flux 
of 1.6x10"^^ W m~^. The total fluxes in the [Arvi] and [Navil] lin es would be 



LutzetaL 



(12000h found a 1.7 times 



9 X 10"^^ W m~^ and 4 x 10~^^ W m~^, respectively, 
higher flux for the [Ar vi] line, suggesting that the distribution of this ion, which possesses 
the lowest ionization energy of the fo ur detected, i s sora ewha t more extended than [S i ix] or 



possesses a large and very faint halo. iPrieto et al. 



(120051) and 



Rodriguez- Ardila et al. 



(120061) 



found that the [Sivi] and [Sivil] lines (with ionization energies of 167 eV and 205 eV, 
respectively) each have halos extending to a radius of 70 pc (1.0"), but they do not provide 
the relative contributions of the core and halo. 

The absence of hydrogen 5—4 recombination line emission is another n oteworthy asp ect 



of the 3-5 fim spectrum in Figure l2j The Br 7 line emission detected by 



Davies et al. 



( )l998l ) from the circum-nuclear ring of H II regions at a radius of 15 -16" has 1-2 o r ders o f 



Tamura et al 



fll99lh . 



magnitude lower surface brightness than that seen at the nucleus by 
Thus the weakness of the recombination line emission in Figure 12] is not likely to be due 
to subtraction of br ight line emis s ion in the offset beam, located 12" north of the nucleus 



just inside the ring. 



Davies et al. 



(I2OO7I ) have mapped the Br 7 line emission in the vicinity 
of the nucleus at high angular resolution and have kindly made their data available to us. 
Our limit for the Br a flux in the central 0.2" xl.O"of <8 xlO"^'' W m"^ is 75 times the 
Br 7 line flux in the same aperture. Assuming Case B conditions this limit corresponds to 
an upper limit on the reddening of 3.5 mag between the wavelengths of the two lines, and, 
for a A~^'^ dependence for the extinction, a visual extinction of no more than ~60 mag to 
the Br 7 line-emitting region. 
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4.2. Solid state and molecular absorptions 



The optical depth of the 3.4 Atm 



the values measured by 



Bridger et al. 



feature rep orted here is i n goo d agreement with 



fll994h and 



Mason et al 



but disagrees with valu es about twice as high reported by 



[ 2004) in la. r ger ap ertures 



Imanishi et al 



fll997h and 



Marco fc Brooks! (120031 ) . As the bulk of the cont i nuurn emis sion comes from t he cen tral 



point source, the differences with 



Imanishi et al. 



(119971 ) and Marco k Brooksl mm are 



difficult to account for in terms of positional variations. In the case of the former paper, 



incorrect location of the short wavelength continuum 
but this does not seem to be the explanation for the 



may be responsi 



e for the discrepancy. 



Marco fc Brooksl (120031 ) resuh. Neither 



we, nor any other o: 
feature reported by 



the above observers ha ve detected the nuclear 3.3 /im (PAH) emission 



Marco fc Brooksl (120031 ). It is possible that some of the features that 



they see are spurious, perhaps resulting from incorrect cancellation of strong telluric lines 
present in the 3.25-3.45 /xm interval. 



Ayani et al.l (|l992l ) and iLutz et al.l (|200J) searched for lines of the CO fundamental at 
resolutions of 250 km s~^ and 120 km s~^, respectively. The latter authors set a 3o" upper 
limit of 7% of the continuum for the depths of individual CO lines. In the 4.57-4.77 /im 
portion of the current spectrum, obtained at 150 km s~^ resolution and also devoid of 
CO lines, the signal-to-noise ratio on the continuum is about three times higher. Thus, 
the limits both on individual lines separated from one another (i.e., with widths less than 
^500 km s^^) and on a broad depression due to blended CO lines are significantly more 

mm . 



stringent than those of 



Lutz et al 



Finally, 



Ayani et al. 



(119921 ) reported the detection toward the nucleus of NGC 1068 of 
a narrow absorption feature at 4.69 /im, which they attributed to solid CO. Our spectrum 
shows no such feature, to a limiting absorption depth of a few percent of the continuum. 
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5. Discussion 

5.1. Excitation mechanisms for coronal line emission 

One of the most important questions concerning coronal lines in AGN is the excitation 
mechanism that drives their emission. Here we discuss whether photoionization or shock 
ionization can account for the locations and strengths of the coronal line emission in 
NGC 1068. A summary is provided at the end of this section. 



5.1.1. Photoionization 



Ferguson et al.l (119971 ) carried out a number of illustrative photoionization simulations 
in an effort to identify the optimal conditions and locations in which the coronal lines 
form around AGN. Assuming ultraviolet radiation from the central engine as the only 
excitation mechanism and plane-parallel, constant density slabs of gas, they determined the 
distances from the ionizing source in which the lines are emitted as a function of density. 
They provide line equivalent widths in the density- distance plane, indicating where the 
bulk of the emission occurs for each line, and allowing rough comparisons between observed 
and predicted emission line flux ratios. The [Alvi], [Arvi] and [Siix] lines observed here 
are included in their calculations. Their flgures assume an ionizing ultraviolet continuum 
simila r to that o f a typ ical Seyfert galaxy with L\r,n 



L 



1/2 



Pier et al. 



(|l99J) and 



Bland-Hawthorn et al. 



]^g43.5 gj^gg 1^ distances scale as 



(119971 ) estimate ionizing luminosities 



1-2 orders of magnitude higher in NGC 1068. 



In the discussion below we assume an ionizing luminosity of 10^^ ergs ^, which 

implies lo cations and dim e nsion s of the emission regions are flve times larger than those 
I 1 I 1 

shown by iFerguson et al.l (jl997i) . if all other parameters remain constant. The model 



then predicts that the distance and density for peak [Silx] line emission cover 1.5-100 pc 
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and lO^-lO^'^^ cm~^, respectively (with the largest distances corresponding to the lowest 
densities). The observed distance of 20 pc (0.3") from the nucleus to the peak line emission 
corresponds to n^lO^ cm~'^). If the actual peak is closer to the AGN than observed but 
obscured by dust, then higher density gas would be required. 



Ferguson et al. 



(119971 ) models predict 



For the [Alvi] and [Arvi] lines the adjusted 
distances to the peak of roughly 15-400 pc for the same density range. Although the 
spatial distribution of the [Arvi] line flux is less precise than the other lines, it appears to 
be similar to that of the other two lines, with a peak at the same location. Clearly the 
observation that the offset from the AGN to the peaks of these lines of with widely varying 
excitations are the same is incompatible with an ionized medium of uniform density, unless 
the extinction is highly non-uniform across the nucleus. The current observations can 
be roughly reconciled with a slab of density several times 10^ cm~'^ if there is increasing 
extinction of the coronal line emission closer to AGN, of sufficient strength to mask the 
true peak of the [Si ix] line emission. If the coronal line emission region is symmetrically 
distributed to the north and south of the AGN, the extinction would need increase further 
to the south of the AGN in order to block the line emission there from view. Evidence for 
higher extinction to the south is in fact widespread (e.g. 



Mason et al. 


2006; 


Das et al. 


2006 



and references therein). A rough lower limit to the additional extinction required to the 
south based on the assumption of a symmetric distribution of coronal line emission is two 
magnitudes at 3.93 fim, which corresponds to approximately 70 visual magnitudes. If the 
highly non-uniform extinction postulated above is not present, only a density distribution 
of ionized gas strongly peaked 20 pc north of the AGN could result in photoionization 
producing the observed maximum line emission from all of these ions at the same location. 



The calculations of .Ferguson et al 



(119971 ) also indicate that for these ions the regions 



of significant coronal line emission can extend to much larger distances (~100 pc) if gas 
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densities decrease gradually with distance from the AGN. The near agreement of our 
estimates of the total flux from the [Siix] and [Arvi] lines with measurements made by 
others in much larger apertures, imply that the ionized gas density must drop steeply at 
distance greater than ~50 pc from the AGN. 



5.1.2. Shock ionization 



Evidence supporting shock ionization in the narrow line regions of AGNs abounds 



(Marconi et al. 


1996; 


Com 


tini et al. 


2002; 


Prieto et al. 


2005 


In the case of NGC 1068, 


Axon et al. 


( 


1998 


) found that 1.5 



Rodriguez- Ardila et al. 



20061 ). 



enhanced at the location of the radio jet. They interpreted this as evidence of interaction 
between the jet and the surrounding medium. The radio jet proceeds outward from the 



nucleus for an angular dis tance of several arc-se c onds, first nearly due north and t 



len at 



posi tion angle 35 deg rees (iGallimore et al. 



and 



Das et al. 



1996, 



2004 ). More recently. 



Das et al. 



fl2006h 



(120071 ) have studied the influence of the radio jet on the kinematics of 
[O III] AA4959,5007. They observed an increase in the radial velocity roughly proportional 
to distance from the nucleus followed by a linear decrease beyond ~100 pc (1.4") and 
concluded that the magnitudes of these changes were so large that the neither can the jet be 
the sole driving force nor can gravity be the sole decelerating mechanism (frictional forces 
must be significant). However, in the case of that line the deceleration occurs well outside 
of the region where we observe the higher excitation 3-5 /im coronal lines. 

The line emission r eported here origin ates i n ionized regions much closer to the nucleus 



than those observed by 



Axon et al 



(119981 ) and 



Das et al 



(120071 ). Nevertheless, evidence 



for shock excitation of th e coronal lines is also r eadily apparent when one compares the 



5 GHz im a ge in Fig. 1 of 



Roy et al. 



Gallimore et al. 



(120041 ) with our data. Assuming, as they and 



(119981 ) have inferred, that source SI is coincident with the AGN, and hence with 
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the infrared continuum peak, the peaks of [Siix] and [Alvi] hne emission at 0.3" north 
are coincident with the bright radio knot C, where the jet bends to the northeast. This is 
highly suggestive of an interaction between the jet and gas in the NLR. Also suggestive 
of an interaction at this loc ation is a bright knot of H2 2.12 fim line emission found by 



Miiller Sanchez et al.l (l2009l ) connecting two streamers of line emission, one of which they 



interpret as feeding the innermost few parsecs. 

The velocity profiles of the of [Si ix] line along the spatial direction, shown in Figure HI 
also appear to support a jet-NLR interaction. At 0.2" north, a prominent blue shoulder 
is evident in the line profile. That structure also appears to be present 0.4" north, but it 
weakens significantly at 0.6". The shoulder at 0.2-0.4" could be associated with ionized 
gas accelerated toward us by the interaction of the jet with the ambient gas near knot C. 
The line profile variations are somewhat remin iscent of the behavior of the [Fevil] A6087 
line reported by iRodrfguez-Ardila et al.l (120061 ) . but the latter is at a much larger distance 
north of the AGN and, also unlike the infrared lines, far from where that line has its peak 
intensity. 

If shocks are the dominant excitation mechanism for these coronal lines at radio knot 
C, then it seems likely that some of the remaining [Silx] and [Alvi] line flux outside of 
that knot originates in other locali zed regions and i s asso ciated with other radio knots. 



The knot labeled "NE" in Fig. 1 of 



Gallimore et al. 



(I2OOJ), located 0.4" to the NE of knot 



'C" and not included in our slit, is a possible candidate for some of this mi ssing line flux 



iowever, judging from the large extent over which [Sivil] was observed by 



Prieto et al. 



2OO5I ) . we do not think that shocks totally dominate over photoionization, and expect that 



some contribution of the latter is present, perhaps even at radio knot C. 



5.1.3. Summary 



The infrared coronal line emission peaks 20 pc to the north of the AGN for ions of 
a wide range of excitation, and extends to 50 pc in the highly excited [Siix] line and 
presumably at least as far in the other lines. To the south there is no clear detection of 
any of the 3-5 /xm lines. Increased extinction to the south can account for the north-south 
asymmetry, regardless of the line excitation mechanism. However, photoionization by the 
AGN producing peaks at equal northern distances from the AGN for such a large range 
of excitations requires a highly asymmetric gas distribution and/or a highly asymmetric 
distribution of obscuring dust to the north. This explanation appears contrived. Evidence 
for collisional excitation of these lines is more compelling; it consists of the spatial 
coincidence between radio knot C (where the jet changes direction), a prominent knot of 
H2 line emission, and the emission peaks of all of the 3-5 /im coronal lines. The presence at 
the same location of a localized prominent high velocity blue-shifted shoulder on the [Si ix] 
line provides further support. Higher angular resolution studies of the coronal line emission 
would allow a more detailed comparison with the morphology of the radio emission. 



5.2. The 3.4 /im absorption feature and the non-detection of CO 



In his subsection we discuss that portion of the interstellar medium producing 
the hydrocarbon and silicate absorption features at 3.4 /xm and 9.7 /xm, but little or no 
absorption by gaseous CO at 4.5-5.0 /xm. O ur discussion makes use of the silicate absorption 



feature as reported by 



Mason et al. 



(120061 ): its similar behavior across the nucleus as the 
3.4 /zm feature provides a strong clue as to the nature of the interstellar medium. The 
subsection on CO examines a number of possible explanations for its absence, and results in 



a similar conclusion regarding the ISM. A summary is provided at the end of this section. 
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It is important to note that in the direction along the (NS) sht the continuum source 
was partially resolved at UKIRT and therefore we are assured that it has a characteristic 
dimension of a few tens of parsecs. Thus the bulk of the material producing the 3.4 fim 
absorption feature and not producing CO absorption lies 10 pc or more from the central 
engine. Dust and gas in the central few parsecs can only be a minor contributor to the 
present spectrum and our arguments as to the nature of the interstellar medium therefore 
apply to the more extended region and not to the "core" or the innermost region of the 
torus observed via interferometry. 



5.2.1. The 3.4 /xm band 



The spectral profile of the 3.4 fim feature is well characterized in the current data. 
Overall the profile is smooth, showing only the two principal sub-features seen in Galactic 
sources, those at 3.420 fim and 3.485 fim. The profile in NGC 1068 contains no evidence 
of a distinct 3.385 /im sub-feature, which is s ometimes present in Galactic sources an d 



sometimes blended with the 3.420 iim feature (jPendleton et al.l 



1994 : 



Chiar et al. 



2OO2I). 



The data also reveal that the depth of the 3.4 /im feature has a similar north-south 
variation across the nucleus as the 9.7 /im silicate absorption (Fig. [5]). Both features are at 
maximum depth on and just to the south of the nucleus. The similarity is also quantitative; 
the fractional increase in the depth of each feature at its maximum compared to the depth 
at adjacent locations is ~40 percent. This indicates that the carriers of the features are 
significantly mixed. The compactness and nuclear location of the region of enhanced 3.4 fim 
and silicate extinction suggests that the extra absorption to the south occurs quite close to 
the nucleus. In the silicate feature the weaker off-nuclear absorption extends over 2" to the 
north and south. Because the 3 fim continuum source is much more compact, the depth of 
the 3.4 nm feature could not be measured more than 0.4" distant from the nucleus. The 
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lower level and more extended silicate feature largely arises in a circumnuclear region of 
much larger dimensions than the torus-like structures seen with mid-IR interferometry. It 
may be associated with material in the disk of the host galaxy, which crosses the line of 



sight s outh of the nuc . 



sharp (jPackham et al 



eus. The distinction between these two regions is not necessarily 



20071). 



The ratio of the optical depths of the 3.4 fim and silicate features at the infrared 
continuum peak of NGC 1068 is 0.19 ± 0.02. In the Galaxy the silicate feature is found 
in both diffuse and dense interstellar clouds, whereas the hydrocarbon feature is found 
only in diffuse clouds. In diffuse clouds far from the Galactic center the ratio of optical 



depths of the t wo features can 



)e de rived from their opti cal depths relative to A 



19841 ) and Av/t3A2 ~ 250 flPendleton et al 



V, 



19941 ). yielding 



Av/rg.7 «i 18.5 fi Roche fc Aitkenl 
'T3a/'T9.7 = 0.07 ± 0.01, a value roughly forty percent of that measured in NGC 1068. 



Toward the Galactic cent er, where /. 



sources observed by both 



Chiar et al 



30 ma g), the value of 



(12002) and 



d/rgy averaged over 



Roche fc Aitkenl (jl985|) is 0.064 ± 0.017 



the same as the above value to within the uncertainties. The ratio in just the diffus e gas 



on the Galactic center s iqhtline may be estimat ed as follows. 



156 ±16 and 



Pendleton et al. 



mm find 



Roche fc Aitkenl (119851 ) find rg.y =3.6 ± 0.3. However, one-third 



of the visual extinction to the Gal actic center is beliey ed to occur in dense clouds, lo cated 



mostly in intervening spiral arms flWhittet et al 



19971). Following 



Chiar et al. 



(120071 ). for 



values of Ay<12 mag the ratio of visual extinction to silicate optical depth is the same in 
dense clouds as in diffuse clouds, and thus we associate those 10 mag of dense cloud visual 
extinction with a sihcate optical depth of 0.54. This gives T^A/Tg.i ~ 0.23/3.06 = 0.075 for 
Galactic center diffuse clouds, again far less than the ratio toward the nucleus of NGC 1068. 
Even if half of the silicate optical depth toward the Galactic center arises in the dense clouds 
in the Galactic disk, the above ratio would only be 0.15, less than the value in NGC 1068. 
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Clearly then, the 3.4 /xm absorption in NGC 1068 is unusually strong relative to the 
silicate absorption. This cannot be explained away by placing the 3.4 fim absorber in a 
more remote location than the silicate feature, because even in Galactic diffuse clouds the 
optical depth of 0.08 observed in NGC 1068, corresponds to a silicate optical depth ^1, 
significantly larger than is observed. The high ratio could be the result of the silicate 
absorption being significantly filled in, as might arise if thermal emission close to the AGN 
were absorbed in the cooler outer portions of the source. Because such filling in would be 
much less for the 3.4 fim feature, the similar variation in optical depths of the features 
across the n ucleus ten d s to a rgue that this is not the entire explanation, however. As 



suggested by llmanishil (j2000l ). the high ratio could be explained, at least in part, by the 
typical continuum surface at 3.4 fim lying interior to the typical surface at 9.7 /im such 
that the added pathlength, unobservable at 9.7 /im, contains a significant fraction of the 
observed column density of the 3.4 fim absorber. The result also could be a consequence 
of either unusual elemental abundances or the chemistry of the dust in the vicinity of the 
AGN. 



5.2.2. The non- detection of CO 
The CO fundamental band has been searched for previously witho ut success in 



NGC 1068 at spectral resolutions ra nging from 20 km s ^ to 250 km s ^ (lAyani et al. 



Lutz et al. 


2004; 


Mason et al. 


2006) 



1992 



20061 ). The current data, obtained at 150 km s~^ resolution, 
significantly tightens the limit on absorption by CO in the warm and presumably rapidly 
moving gas just outside of the 3-5 /im continuum-emitting dust, whose FWHM along the 
ionization cones is ~20 pc. In its lack of de tectable CO in th is band NGC 1068 is similar 



to the five other Type 2 Seyferts studied by 



Lutz et al 



(120041 ) ■ although the limits on CO 



absorption in those galaxies are considerably less strict. 



- 31 - 



We consider three possible explanations for the lack of strong CO absorption lines 
toward the AGN in NGC 1068: (1) the molecules are distributed among many rotational 
states; (2) the molecules reside in a clumpy obscuring medium; and (3) there is a paucity of 
"dense cloud" molecular gas in front of the AGN. To explain the non-detection of transitions 
of CO, OH and H2CO lines from low rotational levels towards the highly obscured nucleus 
of Cygnus A at radio wavelengths, it has been pointed out that at the temperatures of 
the nuclea r clouds the lowest rotationa l levels may contain only a small fraction of t he 



molecules ( iBarvainis &: Antonucci 



1994 



Maloney et al. 



1994 



Conway fc Blanco 



19951 ) ■ For 



example, in LTE at 300 K roughly twenty rotational levels of CO have populations within 
a factor of ten of that of the most populated level, whereas at 25 K only the six lowest 
levels are so highly populated. It also has been suggested by these authors that under some 
circumstances a strong and compact nuclear radio source could redistribute the rotational 
level populations in the surrounding gas, further weakening the strengths of lines from 



the low rotational levels. 



Impellizzeri et al 



(I2OO6I ) find this explanation unlikely, however, 
because the 19% of their sources in which high-excitation OH lines are detected also show 
absorption from the ground state transition. 

In the case of NGC 1068, whose 3-5 /im continuum-emitting dust, and presumably 
the gas associated with this dust, are at temperatures of several hundred Kelvins, the LTE 
fractional populations in the low lying levels normally observed at radio and millimeter 
wavelengths would be much less than in galactic dense clouds. This cannot be the 
explanation for the non-detection of the CO fundamental band lines, however, as the 
M band spectrum in Fig. [2] encompasses lines with a wide range of lower state energies 
(corresponding to J = 1-30 in the case of the P branch). 

The second possibility is that, although there is a large amount of molecular material 
along the line of sight to the nuclear continuum source, it is concentrated into clumps that 
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cover only a small fraction of the source. Each individual clump produces deep CO lines 
but as most of the continuum source is unobscured, the net effect is much weaker CO 
absorption. Such small cloud volume filling factors have been invoked in order to explain 
the mid-infrared observational properties of Seyfert 1 and Seyfert 2 galaxies in the unified 
model (in particular the lack of sili cate emission in the former and the broad spectral 



range of the far-infrared continuum; iNenkova et al. 



20021 ). Small filling fa ctors for toro ids 



surrounding AGNs are also predicted by the hydrodynamic simulations of 



Wadal fl2007h 



If the detectability of CO absorption is related to the cloud filling factor, then for 
similar total extinction CO should be detected in objects with smooth obscuration (i.e., 
covering the entire nuclear continuum source), but not in those with clumpy obscuration. 



Sirocky et al. 



(120081 ) propose that clumpy and smooth obscuration in AGN can be separated 
on the basis of the strength of the 10 and 18 fim silicate features. Strong CO absorption, in 
gas that is warm and presumably close to the nuclear so urce, has been detected in severa. 



AGN - hosting ultralurninous infrared galaxie s (ULIRGs) fISpoon et al 



2006 



2004 



Geballe et al. 



Shirahata et al. 



2007 



Sani et al. 



20081 ). whose silicate features would classify their 



dust distribution as smooth, whereas in the single Seyfert 2 ULIRG whos e spectrum wou 



20081). 



classify it as clumpy, only a weak CO band has been tentatively detected ( iSani et al.l 
This may suggest that the strength of the molecular absorption in AGN is indeed governed 
by the dumpiness of the obscuring medium, but a stringent test of this hypothesis will 
require a larger sample of objects spanning a wide range of dust geometries. 

Nevertheless, the presence in NGC 1068 of the 3.4 fim feature, whose spatial behavior 
suggests that much of it arises in the same gas as much of the silicate feature, close to the 
AGN, raises immediate complications with the clumpy model. That model has the dust 
and gas local to the AGN residing in dense clouds with a small filling factor, so that the 
silicate absorption, which is strong in each clump, is both diluted and filled in by emission 
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( iNenkova et al.ll2002l ). The CO absorption lines would also be diluted by this arrangement, 
apparently to undetectability in the case of NGC 1068. However, the absence of the 3.4 fim 
band in Galactic dense molecular clouds suggests that in NGC 1068 the carrier of that band 
would not reside in such clumps, and thus that a second, more diffuse medium is needed 
for it (ignoring, for the moment the evidence that the carriers of the two dust features are 
mixed). The diffuse material could, for example, be a low-density interclump medium. 
However, the 3.4 fim band in Galactic diffuse clouds is invariably accompanied by a silicate 
absorption feature. As pointed out earlier, the previously estimated value of t^a/ 'T'9.7 in the 
diffuse interstellar medium toward the Galactic center, when applied to NGC 1068, implies 
that most or all of silicate absorption in NGC 1068 must arise in that material carrying the 
3.4 yum band. 

The by now almost inevitable solution to these problems is the third explanation: 
at 3-13 yum essentially all of the absorbing material in the line of sight has the chemical 
characteristics of Galactic diffuse clouds. In this explanation clumpy absorbing material 
is not needed to account for the lack of detectable CO, because in diffuse clouds CO is 
mostly dissociated, comprising only one percent of the carbon. Diffuse clouds do have 



large molecular components, with H2 accounting for typicall y half of the hydrogen ; 



Miiller Sanchez et al. 



;hus H2 



(120091) 



should be available to provide the line emission observed by 
The dust grains in diffuse clouds are comprised of a refractory component responsible for 
the 3.4 /im feature and a silicate component. Gas densities in Galactic diffuse clouds are 
10-1000 cm~^. Densities in the inner few tens of parsecs of NGC 1068 might be greater 
than this. However, at distances of only a few tens of parsecs from a luminous AGN it may 
well be possible that somewhat denser gas also has the above characteristics. 

As discussed earlier, in the Galaxy fairly tight relationships have been found between 



the visual extinction in the diffuse ISM and the depth of the 3.4 /xm feature (jPendleton et al. 
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19941 ) for sightlines to the Galactic center and to local diffuse clouds. Here we employ the 
former relationship, Av/t^a = 150. Using it we estimate a typical visual extinction of 
12 mag to the 3.4 /im continuum-emitting region, implying for a normal gas-to-dust ratio 
that A^(H) = 2.3 x 10^^ cm~^. Then, assuming that, as in Galactic diffuse clouds, one 
percent of the carbon is in CO, and also (conservatively) that the typical continuum surface 
at 3.4 /xm and 4.7 /im are at the same depth, A^(CO) = 6 x 10^^ cm~^, which is consistent 
with the upper limit of 1 x 10^'' cm~^ set by these observations. The above column density 
of hydrogen is much less th an the column de nsity that attenuates our view of the central 
engine, iV(H) > 10^^ cm~^ ( iMatt et al.ll2004l ). but much of that gas must be interior to the 



dust that is emitting the observed 3-5 /xm continuum. 



5. 2. 3. Summary 

We have found that the 3.4 /xm hydrocarbon and 9.7 /xm silicate bands have spatial 
variations that are similar across the central 0.2 x 0.8 arcsecond of NGC 1068 and that 
only a small fraction of interstellar carbon in this region can be in the form of CO. The 
first finding indicates that the carriers of the 3.4 /xm and 9.7 /im bands are largely mixed 
and that significant fractions of each are located in material very close to the nucleus. In 
the Galaxy both features are present in diffuse interstellar clouds, but the 3.4 /im band is 
only observed there. The fraction of CO in diffuse clouds is small and our observed upper 
limit in NGC 1068 is consistent with such a location. Measurable CO absorption would be 
present if even a modest fraction of the silicate absorption arose in fully molecular material. 
The straightforward conclusion is that the 3-13 absorption spectrum of NGC 1068 is 
produced in material chemically similar to the Galactic diffuse interstellar medium. 

We reiterate the significance of the lack of detection of CO and its incompatibility with 
models employing clumps of dense molecular material, even if such material near an AGN is 
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capable of producing the 3.4 /im feature. If the filhng factor of such fully molecular clumps 
were low enough to "dilute" the CO band to non-detection, as observed, similar dilutions 
would occur for the silicate and 3.4 fim absorptions, contradicting the observations. 



The adaptive optics measurements of 



Gratadour et al. 



(120061 ) show that nuclear thermal 



infrared continuum arises in part from numerous clumps covering a few tens of parsecs. If 
the source of the 3-13 fim continuum is dominated by such clumps, then absorbing material 
could either be located on surfaces of the clumps or be more continuously distributed in 
front of them. However, if a significant portion of the continuum is emitted by an extended 
component, the present observations constrain the degree of dumpiness in the absorbing 
material. In that case a low, e.g. <0.1, filling factor for the absorbing clumps probably 
would be inconsistent with the observed strength of the 3.4 fim feature, as it would imply, 
e.g., r(3.4 fim) > 0.8 in the clumps. The maximum observed optical dep th of the feature is 



0.5, in the deeply buried ULIRG IRAS 08572+3915 flGeballe et al 



20061 ) 



6. Conclusions 

The observations reported here bear largely on two important issues: (1) the ionization 
mechanism for the high excitation coronal line emission close to the AGN; and (2) the 
nature of the dusty environment absorbing the infrared continuum of the AGN. Regarding 
the ionized gas, the evidence appears to both favor collisional excitation of the bulk of the 
coronal line emission and, in view of the spatial coincidences of the various lines, argue 
against photoionization. Regarding the absorbing material, the simplest explanation is that 
it has much more of the characteristics of diffuse clouds than dense clouds (clumped or 
unclumped). The strength of the hydrocarbon feature argues against it being distributed 
in front of the sources of the continuum with a very low filling factor relative to them. 
Because elsewhere it has always been found that silicate dust is always present where a 
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3.4 /iin absorption is observed, one can conclude that if the latter is not so distributed then 
neither is the former. 

Finally, it should be noted that in NGC 1068 the extended source of the 3-5 /im 
continuum radiation has a characteristic north-south dimension of ~20 pc, indicating that 
the bulk of the continuum emission occurs at distances of ~10 pc or more from the central 
engine. This distance is much greater than that estimated under the assumption that the 
temperature of the dust is determined principally by the luminosity of the A GN and its 



distance. Using an AGN luminosity of 1.5 x 10 (jWeinberger et al. 



19991 ). dust grains 



Gratadour et al. 



or temperatures 



(120061 ) has suggested 



heated by isotropic radiation from the AGN to t he observed 3-5 Atm col 
of ~500 K would lie only ~1 pc from the AGN. 
that transient heating of small grains could be responsible for the high temperatures. We 
speculate that at least in the north-south direction mechanical heating of the dust and gas. 



perhaps from the AGN jets, could also be important. 



We thank the staff of the United Kingdom Infrared Telescope, which is operated by 
the Joint Astronomy Centre on behalf of the Science and Technology Facilities Council 
of the U.K. TRG's and REM's research is supported by the Gemini Observatory, which 
is operated by the Association of Universities for Research in Astronomy, Inc., under a 
cooperative agreement with the NSF on behalf of the Gemini partnership: the National 
Science Foundation (United States), the Science and Technology Facilities Council (United 
Kingdom), the National Research Council (Canada), CONICYT (Chile), the Australian 
Research Council (Australia), Ministerio da Ciencia e Tecnologia (Brazil) and SECYT 
(Argentina). ARA's research was partially supported by the Brazilian funding agency 
CNPq (311476/2006-6). We are grateful to N. A. Levenson for useful discussions and to the 
anonymous referee for numerous helpful comments. 



-37- 
REFERENCES 

Axon, D. J., Marconi, A., Capetti, A., Maccetto, F. D., Schreier, E., & Robinson, A. 1998, 
ApJ, 496, L75 

Ayani, K., Nakagawa, T., & Nagata, T. 1992, in "Testing the AGN paradigm; Proceedings 
of the 2nd Annual Topical Astrophysics Conference," ed. S. Holt, S. Neff, & C. M. 
Ury, AIP Conference Proceedings, 254, 528. 

Barvainis, R. & Antonucci, R. 1994, AJ, 107, 1291 

Bian, W. & Gu, Q. 2007, ApJ, 657, 159 

Bland-Hawthorn, J., Lumsden, S. L., Voit, G. M., Cecil, G. N., & Weisheit, J. C. 1997, 
ApSS, 248, 177 

Bohlin, R. C, Savage, B. D., Drake, J. F. 1978, ApJ, 224, 132 

Bridger, A., Wright, G. S., & Geballe, T. R. 1994, in "Infrared Spectroscopy With Arrays: 
The Next Generation," ed. I.S. McLean, Kluwer (Dordrecht), 537-538 

Capetti, A., Macchetto, F. D., & Lattanzi, M. G. 1997, ApSS, 248, 245 

Chiar, J. E., Adamson, A. J., Pendleton, Y. J., Whittet, D. C. B., Caldwell, D. A., & Gibb, 
E. L. 2002, ApJ, 570, 198 

Chiar, J. E., et al. 2007, ApJ, 666, L73 

Contini, M., Viegas, S.M., & Prieto, M. A. 2002, A&A, 386, 399 
Conway, J. E. & Blanco, P. R. 1995, ApJ, 449, L131 

Dartois, E., Marco, O., Munoz-Caro, G. M., Brooks, K., Deboffle, D, d'Hendecourt, L. 
2004, A&A, 423, 549 



-38- 

Das, v., Crenshaw, D. M., Kraemer, S. B., & Deo, R. P. 2006, ApJ, 132, 620 
Das, v., Crenshaw, D. M., & Kraemer, S. B. 2007, ApJ, 656, 699 
Davies, R. I., Sugai, H., & Ward, M. J. 1998, MNRAS, 300, 388 

Davies, R. I., Mueller Sanchez, F., Genzel, R., Tacconi, L. J., Hicks, E. K. S., Friedrich, S., 
& Sternberg, A. 2007, ApJ, 671, 1388 

Ferguson, J. W., Korista, K. T., Baldwin, J. A., & Ferland, G. J. 1997, ApJ, 487, 122 

Gallimore, J. F., Baum, S. A., O'Dea, C. R, & Pedlar, A. 1996, ApJ, 458, L36 

Gallimore, J. F., Baum, S. A., & O'Dea, C. 2004, ApJ, 613, 794 

Geballe, T. R., Goto, M., Usuda, T., Oka, T., & McCall, B. J. 2006, ApJ, 644, 907 

Geballe, T. R., McCall, B. J., Hinkle, K., & Oka, T. 1999, ApJ, 510, 251 

Gebhardt, K., et al. 2000, ApJ, 539, L13 

Gratadour, D., Rouan, D., Mugnier, L. M., Fusco, T., Clenet, Y., Gendron, E., & Lacombe, 
F. 2006, A&A, 446, 813 

Honig, S. F. & Beckert, T. 2007, MNRAS, 380, 1172 

Honig, S. F., Prieto, M. A., & Beckert, T. 2008, A&A, 485, 33 

Imanishi, M., Terada, H., Sugiyama, K., Motohara, K., Goto, M., & Maihara, T. 1997, 
PASJ, 49, 69 

Imanishi. M. 2000, MNRAS, 319, 331 

Impellizzeri, V., Roy, A. L., & Henkel, C. 2006, Proc. of Sci, 8th Eur. VLBI Network Symp., 
ed. A. Marecki, 35 



-39- 

Jaffe, W., et al. 2004, Nature, 429, 6987 

Krolik, J. H. & Begelman, M. C. 1986, ApJ, 308, L55 

Krolik, J. H. & Begelman, M. C. 1988, ApJ, 329, 702 

Lee, H.-H., Bettens, R. P. A., & Herbst, E. 1996, A&AS, 119, 111 

Lutz, D., Sturm, E., Genzel, R., Moorwood, A. F. M., Alexander, T., Netzer, H., & 
Sternberg, A. 2000, ApJ, 536, 697 

Lutz, D., Sturm, E., Genzel, R., Spoon, H. W. W., & Stacey, G. J. 2004, A&A, 426, L5 

Macchetto, P., Capetti, A., Sparks, W. B., Axon, D. J., & Boksenberg, A. 1994, ApJ, 435, 
L15 

Maloney, P. R., Begelman, M. C, & Rees, M. J. 1994, ApJ, 432, 606 
Marco, O. & AUoin, D. 2000, A&A, 353, 465 
Marco, O., & Brooks, K. J. 2003, A&A, 398, 101 

Marconi, A., van der Werf, P. P., Moorwood, A. F. M., & Oliva, E. 1996, A&A, 315, 335 

Mason, R. E., Wright, G., Pendleton, Y., & Adamson, A. 2004, ApJ, 613, 770 

Mason, R. E., Geballe, T. R., Packham, C., Levenson, N. A., Elitzur, M., Fisher, R. S., & 
Perlman, E. 2006, ApJ, 640, 612 

Mason, R. E., Wright, G. S., Adamson, A., & Pendleton, Y. 207, ApJ, 656, 798 

Matt, G., Bianchi, S., Guainazzi, M., & Molendi, S. 2004, A&A, 414, 155 

Miiller Sanchez, F., Davies, R. I., Genzel, R., Tacconi, L. J., Eisenhauer, F., Hicks, E. K. 
S., Friedrich, S., & Sternberg, A. 2009, ApJ, 691, 749 



-40- 

Muxlow, T. W. B., Pedlar, A., HoUoway, A. J., Gallimore, J. F., & Antonucci, R. R. J. 
1996, MNRAS, 278, 854 

Nenkova, M., Ivezic, Z, & Elitzur, M. 2002, ApJ, 570, L9 

Packham, C, Radomski, J. T., Roche, P. F., Aitken, D. K., Perlman, E., Alonso-Herrero, 
A., Colina, L., & Telesco, C. M. 2005, ApJ, 618, L17 

Packham, C, et al. 2007, ApJ, 661, L29 

Pendleton, Y. J., Sandford, S. A., AUamandola, L. J., Tielens, A. G. G. M., & Sellgren, K. 
1994, ApJ, 437, 683 

Pier, E/ A/, Antonucci, R., Hurt, T., Kriss, G., & Krolik, J. 1994, ApJ, 428, 124 

Poncelet, A. & Perrin, G., & Sol, H. 2006, A&A, 450, 483 

Prieto, M. A., Marco, O., & Gallimore, J. 2005, MNRAS, 364, L28 

Raban, D., Jaffe, W., Rottgering, H., Meisenheimer, K., & Tristram, K. R. W. 2009, 
MNRAS, 394, 1325 

Ramsay Howat, S., et al. 2004, Proc. SPIE, 5492, 1160 

Rhee, J. H. & Larkin J. E. 2006, ApJ, 640, 625 

Roche, P. F. & Aitken, D. A. 1984, MNRAS, 208, 481 

Roche, P. F. & Aitken, D. K. 1985, MNRAS, 215, 425 

Roche, P. F., Packham, C., Telesco, C. M., Radomski, J. T., Alonso-Herrero, A., Aitken, 
D. A., Colina, L., & Perlman, E. 2006, MNRAS, 367, 1689 

Roche, P. F., Packham, C., Aitken, D. K., & Mason, R. E. 2007, MNRAS, 375, 99 



-41 - 

Rodriguez-Ardila, A., Prieto, M. A., Viegas, FS., & Gruenwald, R. 2006, ApJ, 653, 1098 

Roy, A. L., Colbert, E. J. M., Wilson, A. S., & Ulvestad, J. S. 1998, ApJ, 504, 147 

Sandford, S. A., AUamandola, L. A., Tielens, A. G. G. M., Sellgren, K., Tapia, M., & 
Pendleton, Y. 1991, ApJ, 371, 607 

Sani, E. 2008, ApJ, 675, 96 

Shirahata, M., Nakagawa, T., Goto, M., Usuda, T., Suto, H., & Geballe, T.R. 2006, in 

"The Central Engines of Active Galactic Nuclei," L. C. Ho & J.-M. Wang eds., A. S. 
P. Conf Ser., 373, 505 

Sirocky, M. M., Levenson, N. A., Elitzur, M., Spoon, H. W. W., & Armus, L. 2008, pJ, 678, 
729 

Spoon, H. W. W., et al. 2004, ApJS, 154, 184 

Tamura, M., Kleinmann, S. G., Scoville, N. Z., & Joyce, R. R. 1991, ApJ, 371, 131 
Thompson, R. I., Chary, R.-R., Corbin, M. R., & Epps, H. 2001, ApJ, 558, L97 
Tristram, K. R. W., et al. 2007, A&A, 474, 837 

TuUy, R. B. 1988, Nearby Galaxies Catalog (Cambridge: Cambridge Univ. Press) 
Wada, K. 2007, PAS J, 59, 711 

Weinberger, A. J., Neugebauer, G., & Mathews, K. 1999, AJ, 117, 2748 
Whittet, D. C. B., et al. 1997, ApJ, 490, 729 

Young, S., Packham, C, Mason, R. E., Radomski, J. T., & Telesco, C. M. 2007, MNRAS, 
378, 888 



-42 - 



This manuscript was prepared with the AAS macros v5.2. 



